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Abstract 

Electron and vibrational dynamics of molecules are commonly studied by subjecting them to two 
interactions with a fast actinic pulse that prepares them in a nonstationary state and after a vari¬ 
able delay period T, probing them with a Raman process induced by a combination of a broadband 
and a narrowband pulse. This technique known as femtosecond stimulated Raman spectroscopy 
(FSRS) can effectively probe time resolved vibrational resonances. We show how FSRS signals can 
be modeled and interpreted using the stochastic Liouville equations (SLE) originally developed for 
NMR lineshapes. The SLE provides a convenient simulation protocol that can describe complex 
dynamics due to coupling to collective coordinates at much lower cost that a full dynamical simu¬ 
lation. The origin of the dispersive features which appear when there is no separation of timescales 
between vibrational variations and dephasing is clarified. 
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I. INTRODUCTION 


Stimulated Raman spectroscopy is a common established time-resolved technique for 

(9- 


monitoring vibrational motions |l|-8]. Multidimensional Raman techniques p-luj probe 
the molecular system at multiple time points via a sequence of Raman processes, which 
measure correlations between several coherence periods. Nomac 
in chemical reactions 


12 


13] as well as structural changes [5 


iabatic relaxation dynamics 
, sl, Q] can then be probed 


with high temporal resolution. In a typical UV-(or visible) pump - Raman probe experiment, 
an actinic pump pulse launches a photochemical process in an excited electronic state, which 
is subsequently probed by a delayed Raman pulse sequence. Several variants of spontaneous 
and stimulated Raman probe-techniques which show high temporal and spectral resolution 
have been reported [ 15 I, 16]. In the femtosecond stimulated Raman spectroscopy (FSRS) 
technique the Raman probe sequence consists of a picosecond pulse £2 superimposed with 
a femtosecond laser pulse £3 which stimulates the Raman signal. Starting from the earlier 
work of Yoshizawa and Kurosawa 17] this technique has shown to be a sensitive local probe 
for ultrafast photo-induced processes jl, 6]. Different configurations of the FSRS techniques 
including temporally and spectrally overlapping pulses and resonant Raman processes 18] 
and cascading effects in FSRS [h]] have been calculated [3]. 

Typically in off-resonant FSRS a spectrally resolved pattern of narrow vibrational lines 
(linewidth ~10 cm -1 ) is recorded in short time intervals (20 fs). FSRS is thus considered an 
ideal probe for ultrafast light-induced processes [f, 6] which relates nuclear rearrangements 
to spectral changes. 

Recently we investigated the microscopic origin of the temporal and spectral resolution 
of FSRS [20]. We proposed three levels of theory which form a hierarchy of approximations 


for the simulation of the matter response m based on loop diagrams in the frequency 
domain. These include sum over states and direct propagation of non adiabatic dynamics 
which is usually described by semiclassical models where bath degrees of freedom are treated 
classically. The simplest level of modeling assumes that the vibrational frequency becomes 
time dependent due to e.g structural rearrangement 6, [2QJ In this case the frequency 
trajectory is inserted into the wavefunction description of the signal and the bath degrees 
of freedom are not explicitly included in the Hilbert space description. A higher level of 
theory assumes that some Hamiltonian parameters are fluctuating due to coupling to a 
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bath described by collective coordinates. In that case the signals can be described using 


the stochastic Liouville equations (SLE) that act in the joint system plus bath space 


21 - 


23]. If the bath is harmonic and only modulates the frequencies one can solve the dynamics 


analytically using the cumulant expansion since the fluctuations are Gaussian, and avoid the 
SLE (see Eq. (39) of Ref. j^J). However a much broader class of models with continuous or 
discrete collective coordinates which may be coupled to any Hamiltonian parameter (not just 
the frequencies) can be treated by SLE. Adding the bath complicates the calculation but it is 
still less demanding than the complete microscopic simulation involving all relevant degrees 


of freedom. The SLE, originally developed by Kubo and Anderson in NMR 


2114231] thus 


provide an affordable and practical level of modeling of complex lineshapes. In our earlier 
work we applied the SLE for comparison of the resolution of several Raman techniques for a 


simple two-state kinetics model 


24], Recently we have applied the SLE for identifying the 


spectroscopic signature of the underlying bond splitting mechanism in cyclobutane thymine 
dimer, one of the major lesions in DNA 25]. In the present work we provide more in depth 
general analysis of the SLE in the context of Stimulated Raman Spectroscopy and show how 
to use the SLE to trace the origin of the dispersive lineshapes reported in FSRS experiments. 


II. TRANSIENT ABSORPTION OF A SHAPED PULSE; LINEAR ANALOGUE 
OF FSRS 


We first present a formal expression for the pump-probe signal defined as the change in the 
frequency-dispersed probe intensity . We consider a multilevel quantum system described 
by the Hamiltonian Hq and coupled to an external optical field by the interaction 

H int (t)=£(t)V^ + £*(t)V (1) 

where £{t) and £*(t) are the positive and negative frequency components of the total electric 
field operator E{t) = £{t) + £*(t) respectively. The dipole operator is V = V + W where 
1/1 (V) is the raising (lowering) operator responsible for excitation (de-excitation) between 
the molecular states. The total Hamiltonian is given by 


H T {t) =H 0 + ) 


( 2 ) 


We use superoperator notation 


26-28] that allows to derive compact expressions for 


the signal. With each Hilbert space operator A, we associate two superoperators, denoted 
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as A l (left) and A R (right) defined through their action on Hilbert space operator X as 
A l X = AX, A r X = X A. We further define the linear combinations of these superoperators 
A .|_ = (Al+A r )/2 and A_ = Ar—Ar. A + (A_) superoperator in Liouville space corresponds 
to an anticommutation (commutation) operation in Hilbert space. Using this notation the 
heterodyne detected signal (frequency dispersed transmission of £ p centered around t = t 0 ) 
can be recast as 


S(co,t — to, r) = —Im 
h 


£*(uj) f dte l0j{t to) (TV L (t) exp ( - f dni/ int _(ri))) 

J —oo J —oo 


( 3 ) 


where T represents the set of parameters of the incoming fields. The angular bracket (• • •) 
represents the average with respect to the initially prepared molecular density matrix. We 
define the interaction picture superoperator as A u {t) = exp (iH 0 -(t — tq))A u exp(— iH 0 _(t — To)), 
v = L, R. We also define the retarded Liouville space evolution operator as Git — To) = 
(—i/h) 6{t — t 0 ) exp [— j^Ho-it — r 0 )] . In the frequency domain the propagator is given as 
G(u>) = dte zu>( ' t ~ T0 ' > G(t — To) = \ix>I ~ \H 0 - + ie) -1 where e is an infinitesimal positive 
number used to satisfy the causality condition. Finally T is the time-ordering superoperator 
that orders superoperators in increasing time argument from right to left i.e., 


TA v (ti)B u f{t 2 ) = 9(ti - t 2 )A u (ti)B u f(t 2 ) + 0(t 2 - ti)B u i(t 2 )A„(ti), v, v' = L, R. (4) 


Specific signals are obtained by expanding Eq. (ED to the desired order in the field. 

This paper focuses on the FSRS signal, which is a six-wave mixing process and thus we 
expand Eq. (|3]) to fifth order in H iut _. Before we proceed to FSRS however, we discuss 
a lower order signal which uses the same probe and may provide similar type of physical 
and chemical information. This is a Transient Absorption of the shaped pulse (TASP) with 
a visible-pump and shaped probe (broad plus narrow band) which is a linear analogue of 
the FSRS and is obtained by expanding Eq. (j3J) to third order in H , nt _. The measurement 
is the transient absorption (TA) of a shaped pulse consisting of a narrow plus broadband 
probe (Fig. 1). An actinic pulse £ a , centered at tq, promotes the molecule from the ground 
electronic state \g) to a superposition of vibrational levels of an electronic excited state. A 
spectrally narrow £\ and spectrally broad £ p probe pulse, both centered around t 0 > r 0 , 
then interact with the molecule and can either stimulate emission to lower vibrational state 
or absorption to a higher vibrational state and the frequency dispersed transmission of the 
probe pulse is recorded. The positive frequency component of the electric field is given by 


£{t) = £ a (t - r 0 ) + £i(t- t 0 ) + £ p {t - to). 


( 5 ) 
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FIG. 1. TASP signal: Level scheme (a) and ladder diagrams (Eq. (j6j ) and Eq. flj]) ) (b). 


The entire sequence of events (including the actinic pulse) is a four wave mixing process and 
the signal scale as £i£ p £%. The signal can be viewed as a generalized linear response with 
respect to £\£ p from the non-stationary state prepared by £ a . 

Fig. 1 represents the level scheme and corresponding ladder diagrams that contribute to 
the signal (for diagram rules see ref. j^)- We first expand the exponential in Eq. (EJ) to first 
order in £\ and £,, and to second order in £ a , taking the actinic pulses to be impulsive £ a (r) = 
£ a S(r). The following expressions can be directly read off the diagrams (see Appendix A) 


°TASP 


(u;to,T 0 ) 



-) dt dt 3 e luJ ^~ to) (V L Q{t 3 )V ] R G{t-t 3 -T 0 ))' 
' l J og J 0 


x £*(u)£i(t-t 3 -to) , 


( 6 ) 


n(ii) 

°TASP 




T Im 

h 




dt 


dt 3 (V L g(t 3 )Vlg(t-t 3 -T 0 ))' 


x S*(u)£i(t-t 3 -to) 


(7) 


where (■■■)' = JZ ac Pac{(I\ • • • | ac)). If the system is prepared in a superposition of vibra¬ 
tional states a and c in an electronic excited state, created by two interactions with actinic 
pulse £ a represented by the action of VrV^ on the ground electronic state. A non-stationary 
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vibrational wave-packet in the excited state [2| is then described by \£ a \' 2 Vtf\ / l\gg)) = 
Pac\ ac )) where p ac = \£ a \ 2 PgaPcg- By inverse Fourier transformation of the propagator 
and the electric field, performing time integrals and introducing a delay between preparation 
time To and probe pulse to: T = to — tq Eqs. jBJ) - (Ej) yield 

O • pOO 7 A 

■ W ( - rr ^~ ' ~- iAT S;(uj)£ 1 (uj + A)x ( T ) ASP{i) (-u,u + A) 


QW , 2 r * f°° dA _ 

Stasp^T) -_ -e 


StASp( U i T) — ^ I m 


L<I> 


/ —OO 

dA 


e iAT £*(uj)£i{ll> + A )Xta SP (u)(- uj ’ uj + A ) 


(i) 


2 n 


( 8 ) 

(9) 


where we have introduced the generalized susceptibility Xjytsp(j )( —Iu + A), j = i,ii 

x i T ) Asp i p(-“^ + A) = (v L g(u)v£g(- a))', (io) 

xfi SP(ji) (- W)W + A) = < V L Q(u)VlQ{-A ))'. 


( 11 ) 

Four-wave-mixing such as transient absorption was proposed by P. Champion’s group 29, 3p| 
where it has been called “effective linear response approach”. Here we extend the idea to 
recasting the n + m-wave mixing signal in terms of the effective m-wave mixing. We use 
this approach in the context of six-wave mixing FSRS signal. Eqs. (JHJ) - (EJ) are similar to 
those obtained by Champion, but for TASP - the linear analogue of the FSRS. 

The signal may be interpreted as a generalized linear response to the field £\£ p from 
a nonstationary state p ac prepared by the actinic pulse £ a . Note that Xtasp depends ° n 
two frequencies (rather than one for systems initially at equilibrium). This implies that 
the signal (Eqs. fl8j |9])) is sensitive to the phase of the field. The choice of broadband £,, 


and narrowband £\ is one example of more broadly defined pulse shaping 


31-5- 


Different 


phase shapes can manipulae signals by enhancing or suppressing various spectral features 
and changing the line shapes, etc. In our earlier work we have been studying the pulse 
shaping of linear signals as a possible tool for coherent control ( 35 ] . Our formalism can 
incorporate arbitrary pulse shapes which can and provide many novel control tools for the 
signals. 

Assuming that the £\ pulse is spectrally narrow (picosecond) i.e., £\{u j + A) = £\ 5 (u> + 
A — (jJi) we can perform the integral over u>[ which results in 


s tasp(^ T ) = t I m (—i)e <(w ul)T £;(u)£iX { T ) ASP(i)(-UiUi) 


h 


Srlspi^'i T) — ^ I m 


(i) 


) ei(w U1)T ZpH^TASPwi-^iUi) 


. K h 


( 12 ) 

(13) 
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FIG. 2. Ladder diagrams for FSRS signal Eq. (TT7T) and Eq. (USD. 


The total signal is finally given by S^sp ( u > T) = S P \ SP (uj ; T)+S , ^ sp (c(j; T). In the absense 
of a bath the correlation functions can be expanded in terms of molecular eigenstates which 
gives the signal as 


Stasp( r ) 
s tasp( r) 


2 

h 

2 

h 


Im I] 

a,c,d 

Im 52 

a,c,d 




e <(w Wl)T HdaHcd£p(u)£lPac 

^ad T ) (^ kfi ^ac T A'«r) 

e-^-^T ii ad n dc S*(uj)£ip ac 
^da T ^Tda)(^ ^1 ^ca T ^Tca) 


(14) 


These results will later be compared with the FSRS signal Eq. (EHJ). 


III. OFF-RESONANT FSRS 


FSRS is a powerful technique for studying photophysical and photochemical processes 


in molecules 


36] including vibrational and nonadiabatic electron dynamics . It is a six- 


wave mixing process [24] that, when the actinic action is treated implicitly, is given by a 
generalized Moreover in the electronically off resonant regime, considered here, it can 
be expressed in terms of a generalized linear response y- 1 ^ that resembles TASP. 

The general expression for a four wave-mixing probe of a nonstationary system is given 
in Appendix B. FSRS (Fig. 2) can be described using diagrams similar to Fig. (1) where we 
simply replace the dipole operators by the polarizibility. We use the same level scheme as 
in Fig. (1). The effective radiation-matter interaction Hamiltonian (in the Rotating wave 
approximation) for the electronically off-resonant Raman process induced by pump pulse S\ 
and the probe pulse S p , reads 


= a n £*(t - t 0 ) Siit - t 0 ) + S a (t- r 0 ) V e f + h.c. 


(15) 
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Here the dipole operator V e couples the molecule with actinic pulse £ a , centered at time 
To and creates a non-stationary state in the excited electronic level. a n = a n + o£ n is the 
excited state polarizibility that couples parametrically the pump £\ and the probe fields £ p . 
Both pulses arrive simultaneously and are centered at time t 0 > r 0 . Similar to TASP we set 
the delay between the actinic and Raman pulses as T = t 0 — r 0 . Note that this off-resonant 
Raman process is instantaneous and the system can only spend a very short time in the 
intermediate electronic state. Following Eq. (EJ), the frequency-dispersed probe transmission 
signal is written as 


Sfsrs(w, T) — — Im 


/ dt 

' —oo J —oo 


/OO r\, 

—1 ) fi-to) £* (cy) ( aL ( t)e~ m J - Hint ~ (r)dr ) 

-oo 27r 


(16) 


The picosecond pump pulse £\ has a narrow bandwidth centered at frequency uq. This 
signal can be read off from the diagrams given in Fig. (2) and is given as 


S FSRs(u\toiTo) = - Im 


j poo poo 

(t) / dt dt 3 e luJ ^~ to) (a nL g(t 3 )a nL g(t-t 3 -T 0 ))' 

J — oo J 0 


x £p(uj)£ p (t - f 3 - l 0 )£i(t)£r(t-t3-£o) 


(17) 


SFSRsi^t0,T 0 ) = |lm (~) 


POO PC 

/ dt 

1-00 Jo 


) 

dt 3 _t ° ) (a nL g(t 3 ) a nR g(t-t 3 -r Q ))' 


x £*(u)£p(t -t 3 - t 0 )£i(t)£i(t-t 3 -t 0 ) 


(18) 


Eqs. (R7I) - (IT%j) for FSRS or Eqs. (jnj)-(lZI) for TASP can be used for high level numerical 
simulations including nonadiabatic dynamics. Alternatively we can recast the signal using 
frequency domain Green’s functions. Denoting T = t 0 — r 0 we get the total signal as 
^FSRsi 10 * T) = Sp > SRS ( a;; T) + Spg R3 (u>] T), 

e lAT \£i\ 2 £*(u)£ p (u + A)xpg R g^(—(jj+ijJi,—Ldi+Lj+A) , 

(19) 


S^SRsi^ i^ 1 ) — 


Sfsrs( u ’T) — ( ^) 


dA 

2 n 


,-iAT|^ |2$7* 


11 2 £p (w)£ p (uj + A)Xfsrs(u) ( — <^+^11 — uq+aH-A) 


( 20 ) 


The generalized FSRS susceptibility is now given by 

X^F 3 RS(i)(~ ^+^ 1 ) — uq+cj + A) = (a n x ,£/(id — cq)cKnz,f/(—A)) , (21) 

X-FSRS(ii )( k’+kfii (UiTtuTA) k-h)oqiif/( A)^ • 


( 22 ) 













This signal may be viewed as the linear response to £* £\ with the generalized susceptibility 
formally similar to TASP fEqs. (fTT)l) - (fill) ) except that the dipole operators V are replaced 
by the polarizability a n . The reduced description of the six-wave mixing FSRS in terms of 
an effective four-wave-mixing has been studied previously 37] where the response function 
has been calculated using the multimode Brownian oscillator model. The Green’s functions 
expressions presented here provide more general framework for calculating the signals. In 
particular our approach is very effective if the excited state coherence is long lived on the 


time scale of the Raman probe sequence [7, 38(] It also applies in the presence of a non¬ 


trivial relaxation process such as internal conversion between excited states 


37]. In the 


absense of a bath we expand the signals ( 1121 ) - ( 1201 ) in the molecular eigenstate basis and 
obtain 


^FSRsi^’T)— ^ 

a,c,d 


s 


(ii) 

FSRS 


U-T) = -Im J2 

a,c,d 


J f°° dA e- iAr a da a a ^;(^)4((u + A)|£i| 2 p ac 
J— oo ^ad + A UJ ac A i^/ a c) 

,i, r dA e~ lAT a ad a dc £*(u)£ p (u A A)\£i\ 2 p ac 

J—oo {ui UJ\ UJ da A ^yda)( A UJ ca A i^fca) 


■ (23) 


Comparing FSRS signal (Eq. (1231) ) to the TASP (Eq. (fl4l) ). we note several differences. 
First the dipole moments in TASP are replaced by corresponding excited state polarizabili¬ 
ties in the case of FSRS. When the narrowband pulse £\ is monochromatic the TASP signal 
has a perfect frequency resolution, which governs the Raman resonance u — uj\ ~ uo da . In the 
same time it loses all temporal resolution, since the exponent e -*A-^i) T j n Eq m has no 
dependence on material parameters such as energies and dephasing rates. In order to obtain 
some temporal resolution, £\ must have a finite bandwidth. Generally the temporal and 
spectral resolutions of the TASP signal will be Fourier conjugates, so the perfect time reso¬ 
lution would correspond to poor spectral resolution. In contrast, the FSRS signal contains 
an extra integration over A. Therefore, the ac state in Eq. (1221) does not give a spectral 
signature but rather controls the temporal resolution. Furthermore, once we fix the narrow- 
band £\ with frequency uj\ the Raman resonance u — uoi ~ uj da is now well resolved. Thus 
under the same conditions FSRS has both high spectral and temporal resolution, which are 
not directly Fourier conjugates of each other. In the next section we show how the coupling 
to a dynamical bath can be incorporated to calculate signal using the stochastic Liouville 
equation. 
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IV. SIMULATING FSRS SIGNAL BY THE STOCHASTIC LIOUVILLE EQUA¬ 
TIONS 


The stochastic Liouville equations 
for computing spectral line shapes [21 


SLE), first developed for NMR, is a convenient tool 
22 I - 39-41]. This approach assumes that the quantum 
system of interest is affected by a classical bath, whose stochastic dynamics is described by 
a Markovian master equation. The SLE is an equation of motion for the joint system plus 
bath density matrix p 

^ = Cp(t) = -jJH,p(t)] + Lp(t). (24) 


where L represents the stochastic Markovian dynamics of the bath. Both discrete V-state 
jump and continuous collective coordinates (Fokker Planck equations) are commonly used 
to model the bath. In our model, the system has two vibrational states a and c, and the 
vibrational frequency 0 J ca}S is perturbed by the bath which has N states, s = 1, 2, • • • N. The 
entire density matrix p thus has 4 N components \mSs)) which represent the direct product 
of four Liouville space states zv/)), where v, z/ = a,c, and s represents N bath states. The 
Liouville operator C is diagonal in the vibrational space, and is thus represented by four 
N x N diagonal blocks in bath space. 


\£\ w's,vivas' S,l/v'SJ 1 (25) 

where Ls = —K describes the kinetics given by the rate equation: 

= (26) 

s' 

where Paa it) is the population of the s-th bath state. The solution of Eq. 'ft ill is given by 


= £ U„, exp (27) 

s' 

where U is the transformation matrix, where the eigenvectors are organized as rows. This 
matrix satisfies left-eigen equation ^ U sp K ps > = K^fU ss ’ as the rate matrix K is not 
Hermitian. pa} (0) represents the population of the initial bath state. 

The coherent part Cs = —(i/h) [H $,...], which describes the vibrational dynamics, van¬ 
ishes for the |ao)) and |cc)} blocks, [ Cs\aa,aa — [£-s\cc,cc = 0. The remaining blocks of Cs 
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read 


[■^S']ac,ac i 


uj^ 0 ... 0 \ 

o Ul ] ... o 


\ o o ... Ua ] ) 


(28) 


The two Liouville space Green’s functions (i.e., the solution of Eq. are thus given by 


Qaa,aa(t) = ~~9(t)ex p [L S ]t 


= —jQ(t)U exp [[Lsf-H 


U~\ 


(29) 


Qac,ac(t) = - exp ([L S ] + [C S ]ac,ac)t 


= eX P [^lac^act 


v-\ 


(30) 


where U and V are transformation matrices, which diagonalize the matrices in the exponents. 
The time domain FSRS signal on the Stokes side (w < ojq) is given by 


1,00 dA ( 

2 n ip ' 


S fsrs (u>, T) = S I A_£»£> + A)SP srs (u, T; A), 


(31) 


where Sj? SRS ((d, T; A) can be recast in Liouville space as follows: 


r\ pOO pt 

SpSRsi^’T'jA) =— I dt I dr 3 |^i| 2 |^ a | 2 e _jA(T3_T) 


’ — OO J — OO 


x - T 3 ,T 3 ), 


(32) 


where by using the Green’s functions in Eqs. (1291) and (l30l) . the matter correlation function .F(£i ? £ 2 ) 
is given by 


F{tl, t 2 ) = - ^ ^2 a lc\ V ag\ 2 ((I I Gac,ac(tl)Gaa,aa(t2) |po))s 


h 


) 3 E a -lK 9 | 2 ^(ti)^(^)e-^( tl+2t2) 


a,c 


X (1,1,..., 1)Vexp V 1 U exp 


[L s } dias h 


u 


-1 


/1 \ 

0 

Vo ) 


( 33 ) 
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Here, the initial state is the direct product, 


\po))s = I aa)) 


/ 1 \ 
0 

V o y 


(34) 


and we have traced over the final state ((/1 = (1,1,..., l)Tr where Tr = ((aa | + ((cc\. Vibrational 
dephasing terms have been added; e~ lat is added to Gac,ac and e~ 2lat to G a a,aa • Inhomogeneous 
broadening can be included by convoluting the present results with a spectral distribution. Gaus¬ 
sian frequency fluctuations can be incorporated via the cumulant expansion which solves the SLE 
for a Brownian oscillator bath. We have incorporated this level of theory in an earlier study 
Upon evaluating the time integrals in Eq. ([321) we obtain 


20 ]. 


Sfsrs&’T] A) — 


-2 i 


K 2 


l^ 1 | 2 |^| 2 E««clK s | 2 e iA P(/|W(^-^ 1 )^(-A)|po))5. 


(35) 


a,c 


Following Eq. 4ZZD we introduce the bath population of the state a after interaction with the 
actinic pulse: 

P$(t) = \£a\ 2 \Va 9 \ 2 Gaa,aa(t)\P0)) S . (36) 


Substituting Eq. fl36l) into the signal expression Eq. (l32ll and Eq. (l3Tll gives 


— 2 ? 

Sfsrs(u,T) = \£] 


iVa 2 

I / j lJC ac 
a,c 


^ ) Gac,ac,s(bJ 


, [°° dA 


£ p (cu + A)e lAT p^(-A), (37) 


where p aa (— A) is the Fourier transform of the population of the state a. Q a c,ac(w) is a frequency- 
domain Green’s function and represents the sum over bath states. It follows from Eq. (|37l) 
that the A integration represents a path integral over the bandwidth corresponding to the inverse 
dephasing time scale. This integral is generally a complex number. Therefore, the signal dSZD 
depends on both real and imaginary parts of the coherence Green’s function G a c,ac(u\ and thus 
contains absorptive as well as dispersive spectral features. In the limit of slow fluctuations, one can 
neglect the jump dynamics during the dephasing time. In this case we can replace £ p (cu + A) ~ 
£ p ( a;), the integral over A simply yields p aa (T ) and we obtain the static averaged signal 


Sfsrs(u,T) = EE & FSRS,a( U )p { a S ^T), (38) 

a s 
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with 


^FSRS,a( U ) ~ ^^9 l^p( a; )| 2 l^ll 2 YagcOvMsi" UJl )’ (^9) 

c 

For comparison we give the corresponding TASP signal: 

Stasp(u,T ) = XI W 2 X^«c, a c, s (w)e* (w “ Wl)r /3^(a; - wi). (40) 

a,c s 

Unlike the general FSRS signal (I57|) , the static averaging limit (1581) only contains absorptive line 
shapes since bath dynamics is neglected during the dephasing time. Furthermore the time evolution 
in this case is governed by a snapshot of the populations of the excited states. Eq. ([57|) and (1551) 
are therefore expected to be different at short times and become more similar at longer time. 

Note that since the signal (I57|) is written in terms of Green’s functions expanded in sum over 
states, it can be applied to complex systems with multiple vibrational modes coupled to various 
baths. In the typical chemical reaction, such as isomerization, only few collective coordinates are 
involved. Therefore these degrees of freedom are typically treated explicitly whereas the rest of 
the vibrational and bath degrees of freedom can be treated by a harmonic approximation. In the 
following we use a simple model system to illustrate the power of SLE approach . Using Eq. ([571) 
we performed simulations of the signals for a model system with four vibrational modes and 
ten bath states (N — 10). Two vibrational states (i.e., a and c) are included for each mode, 
there are two vibrational states (i.e., a and c). We use a kinetic model described by rate 
equation ([261) and we assume a linear chain of forward and backward reactions among the 
bath states. 

State 1 ^ State 2 ^ ^ State 10. (41) 

k -1 k—2 k -9 

The rate constants vary linearly along the chain 


(42) 


k, = h + - 1 ), 

k-i = 0.1 hi, 

where k\ > Aq. These are given in Table |U The process slows down along the chain. This 
allows to observe both fast and slow jump modulation regimes. Considering the detailed 
balance relation, the second line in Eq. (1421) indicates a constant energy difference between 
two neighboring bath states, namely s and s + 1 states. The resulting population dynamics, 
obtained from Eq. (j27jh is depicted in Fig. SI ( 42 ] . At T = 20 ps, not only the state 10 
but also several bath states contribute to the signal. As shown in Fig. El in our model, the 
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frequency for a given vibrational mode depends on the bath states (Eq. (J2EJ)), and satisfies 
a linear relation. 


(jj 


0) 


— + ^ca(s — 1 ), 


(43) 


where S ca represents the frequency shift when the bath transits from the s to the s +1 state. 
As shown in Fig. [3] and Table HI we employed two parameter regimes both using two values 
of S ca ; small hi is used for two modes (mode 1 and 2, hereafter), and large S 2 is used for the 
remaining two (mode 3 and 4). The two regimes correspond to different relation between 
the jump rate k and splitting 5 ca . In the first regime (I), modes 1 and 2 rapidly evolve ( 
k\ A hi, fast modulation limit - FML) and the modes 3 and 4 are modulated slowly ( k\ -C 6 2 , 
SML). In contrast, in the second parameter regime (II), all four modes are subject to FML. 
Note that at longer delay times, the jump rate itself slows down (£q > kg) as mentioned 
above. The cola of mode 2 and 3 are set so that their u>ca frequencies show crossing between 
each other (see Fig. [3] for the crossing frequencies). We take £\{t) = £ie~ iuJl ^ t ~ T " > to be 
monochromatic, whereas the Raman probe has a Gaussian envelope with center frequency 
oj p and finite duration ex, £ p (t) = g pe -( t ~ T ) 2 / 2rT2 - lu, p( t - T ) _ The integrations of Eq. (ET71) are 
then performed analytically to get Sfsrs^iT), where prefactor u 2 |£i| 2 |£J 2 |£ p | 2 a; 2 c |Vj, 


ag\ 


is 


set to be 1. Note, that in reality the duration of the £\ is finite. However since it is in the 
picosecond range, whereas all dynamical processes as well as the duration of the £ p are in 
the femtosecond range, the the CW approximation is justified. In more general case, the 
bandwidth of £\ has to be taken into account which will reduce the spectral resolution of 


the signal. As has been shown in our earlier work [2Q] the duration of £ p has to be optimized 
in order to obtain both high temporal and spectral resolution. 

The resulting FSRS signals obtained for regime I are shown in Fig. 01 In mode 3 and 4, 
corresponding to SML, fine structure features of bath states are found even at early delay 
times (around 1 ps, see Fig. 0b). The fine structure is reduced for longer delay times (> 10 
ps), reflecting the population decay shown in Fig. SI [42]. This indicates that the signals 
are in the snapshot limit. The SLE coincides with the static average limit. In modes 1 and 
2, fine structures are found only after 5 ps. At T = 2 fs, dispersive lineshapes are present 
since the dynamics is not negligible during the dephasing (see Eq. (j3HD b The FSRS signals 
for regime II are shown in Fig. [5l At early delay times, the signals have no fine structure 
(Fig. Eb). Only in mode 3 and 4 at long delay times (> 5 ps), there are several weak fine 
structure features due to approaching the SML. These signals are beyond the snapshot limit. 


14 





(a) (b) 

FIG. 3. Model frequency change of four vibrational modes along the sequential reaction, (a) 
Parameter regime I and (b) regime II. In Eq. (j43l) . the of mode 2 and 3 are set so that their 

(s) 

(jJca frequencies show crossing during the transition from bath state 5 to 6. The frequencies at the 
crossing points (dashed lines) are shown. 


TABLE I. Parameters employed in the FSRS simulations. 


Parameter regime 

h (s- 1 ) 

kg (s b ) 

Si (s' 1 ) a 

$2 (S' 1 ) a 

(7 (fs) b 7 a (s b ) c 

I 

1.00 x 10 12 

0.667 x 10 12 

3.76 x 10 12 

7.51 x 10 12 

20.0 1.88 x 10 ] 

II 

1.00 x 10 12 

0.333 x 10 12 

0.939 x 10 12 

3.76 x 10 12 

30.0 1.88 x 10 ] 


a The frequency shift S ca is small 5± for two normal modes, and large S 2 for the rest two modes. 
b The center frequency of the Gaussian Raman probe, cu p , is equal to — 1000.0 cm -1 . 
c The vibrational dephasing time corresponds to linewidth of 10.0 cm -1 , and is used for all four 
vibrational modes. 


Furthermore, as the Raman signal is measured using heterodyne detection which involves 
terms that correspond to interference between signal and background contributions. Unlike 
homodyne detection which yields only absorptive line shapes, these interference terms may 
result in dispersive line shapes. The dispersive line shapes, however do not always show 
up in the heterodyne detected signals and their existence depend on the parameters of the 
system as well as laser pulses as we showed above. Therefore the SLE allows to explain the 
dispersive line shapes observed when the dynamics of the system is fast compared to the 
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dephasing time scale, that it cannot be neglected. Such simple analysis is only possible in 
Liouville space and cannot be done in Hilbert space. 




(a) (b) 

FIG. 4. Variation of simulated FSRS signals with different delay times T up to (a) 15 ps and (b) 
0.95 ps. Parameter regime I was employed. After T = 2 fs, time intervals are (a) 500 fs from 
T = 500 fs up to 10 ps, and 1 ps later and (b) 50 fs from T = 50 fs up to 1 ps. The stick spectra 
on the horizontal bottom (top) axis represent the frequencies of the state 1 (state 10). 


V. DISCUSSION AND CONCLUSIONS 

In summary, we have employed a superoperator diagrammatic technique to calculate 
the linear and third order frequency-dispersed transmission signals. Assuming impulsive 
(femtosecond) actinic pulse we show that both TASP and FSRS signals can be recast in 
terms of a generalized linear response function with respect to pump and probe field from 
a non-stationary state prepared by the actinic pulse. Our expressions can be used in three 
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FIG. 5. Variation of simulated FSRS signals with different delay times T up to (a) 15 ps and (b) 
0.95 ps. Parameter regime II was employed. After T = 2 fs, time intervals are (a) 500 fs from 
T = 500 fs up to 10 ps, and 1 ps later and (b) 50 fs from T = 50 fs up to 1 ps. The stick spectra 
on the horizontal bottom (top) axis represent the frequencies of the state 1 (state 10). 


types of simulation protocols: Eq. (JTJ) and Eq. (1TH1) may be used in direct nonadiabatic 
propagation of the wavefunctions. The sum over states expressions Eq. (TT4l) and Eq. (1231) 
may be used for the interpretation of resonances. Finally the stochastic Liouville equations 
(Eq. rmTTRl) provide an inexpensive level of modeling that is numerically less demanding than 
the direct propagation. This approach provides a qualitative explanation for the dispersive 
features observed in Raman signals. It further covers the fast dynamical behavior which 
cannot be described the snapshot limit, which yields the system response far beyond the 
instantaneous frequency tracking. 

We first considered the linear analogue of the FSRS - a transient absorption of the shaped 
(broad plus narrowband) pulse Stasp(u,T) signal: linear transmission of a single probe 


17 








































following a femtosecond actinic-pump pulse. The entire process (including the actinic pulse) 
is a X (3) which can be interpreted as a generalized yA when starting in a nonstationary 
state. The signal with a second order (three wave mixing) probe is a yA process which can 
be recast as a generalized y (2 h For a third oder - four wave mixing (FWM) probe the yA 
process can be described as a generalized yf'T The generalized n-th order susceptibilities 
X <n ) depend on n + 1 rather than n independent frequency variables due to the lack of time 
translational invariance of the time evolution of the system. 

The present formalism can be generalized to other spectroscopic techniques. The off- 
resonant FSRS can be extended to the resonant case, which generally provides a stronger and 
more complex signal. In particular the six-wave mixing process which generally described 
by will be reduced to a four-wave mixing described by a generalized y® compared 
to a generalized x ^ °f a combined broad plus narrowband pulse in the off-resonant case. 


The Impu 
E, 43, 


sive excitation may include photoexcitation by a short optical or infrared pulse 


44], photoionization 45] that prepares an excited ion and core excitation followed 


by an Auger process 


46] that leaves the molecule in a doubly ionized state. 
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APPENDIX A: DERIVATION OF THE TASP SIGNAL 


We read the signal from the diagrams of Fig. 1. First we expand the exponential in 
Eq. (EJ) to the first order £\ and £ p and second order in £ a which yield for the first diagram 

(i) 


^TASP ! ^0? r o) ~ I m 


r*r 3 


r*T2 


dt 


dr 3 / dr 2 / dr^e 1 ^ 1 to) £*(u) £i(r 3 -t 0 ) 


’ —oo J —oo 


(V L g(t-T 3 )V} i Q(T 3 -T2)V R g(T2-T 1 )Vl j ) £*(t 2 -T 0 ) Eafa-To) 

+ (V L g{t-T Z )Vlg{T Z -T 2 )Vlg(T2-Ti)V R ) £ a {T2-T Q )£* a {Ti-T 0 ) . (44) 
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where two terms indicate two possible ways that actinic pulse can excite the ground state to 
the vibrational state of excited electronic state. For the first (second) term the absorption 
happens first on the ket (bra) side followed by the absorption on the bra (ket) side. Now 
assuming the actinic pulses to be impulsive S a (r) = S a 8(r), we get 


STASP^t o,T 0 )=|lm (-^) 


POO POO 

dt / dU e iuj{t ~ to) 


<’ V L G{h)v' R G{t-t z -T 0 )V R vl) 

+ (14^(t3)^(t-t3-ro)^)|^H^(t-t 3 -to)|4| 2 l, (45) 


where note that the first propagator Q from the right is evaluated at t = 0 and we use 
C?(0) = —i/2h. The action of VrVI or V[ V R on ground electronic state creates vibrational 
wave-packet in the excited state |oc)) which is a non-stationary state. Now performing 
inverse Fourier transformation for the propagator and for the electric field and integrating 
over t and t\ variables we obtain the signal as 

duj 0 f°° dui[ ^ i ^ (w / +Wo )^( Wo ))) , e -h--c,Oto e ^oro 

(46) 




(0 


TASP 


2 r i 

(uj]t 0l T 0 ) = — Im|^(——) 


27r ,/_™ 27r 
x T*(o;)Ti(cc; , 1 )|T a | 2 27r(5(a; - ui 0 - u[) 


Finally integrating out u> 0 variable we obtain 

2 „ r, i x r°° dui'i 


S' 


TASpi^t o- T o) = |!m[(-^) 


2 n 


(VlQMV, JS(w-w;))>'£ p *M£ 1 K)|£„| 2 e- i( “-“i> ( '”-’»> 

(47) 


Note that the above expression depends on the time delay between the actinic and the probe 
pulses which we write as T = to — r 0 . Therefore the signal can be recast as 

2 „ r, i s r°° dA 


STASp(W iT) — I m ( J-) 


2 n 


e 1 £*(uj)£i(u> + A)x { tasp ^(-u,u> + A) 


(48) 


where the generalized susceptibility XrASP(i)(~ 00 > u + A) is given as 

+ A ) = <ue(«)v^s(-A)>'. 

The expression for diagram (ii) of Fig. 1 can be derived similarly. 


(49) 


APPENDIX B: THREE AND FOUR-WAVE-MIXING PROBES 


Here we consider a three wave-mixing (TWM) and four wave-mixing (FWM) signal. After 
preparation of the state p ac we sent pulses £\, £•> for TWM, and additional £ 3 for FWM 
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which are centered around Ti, To, and T 3 relative to the preparation time To. We then detect 
the frequency dispersed transmission of the probe pulse S p centered at t = T with respect to 
t 0 . The signal can be calculated in the second and third order of the field matter interactions 
for TWM and FWM respectively and we obtain 


Stwm{u\ T,Ti,T 2 ) — — Im 

n 


f°° AAA g-iwT+iwiTi+ic^Tb 


i—oo J—oo 


x ^p( w )^i( w i)^2( w 2)XrvrM( uj,uji,uj 2 


where the second order generalized susceptibility is given by 


(50) 


(51) 


Similarly the FWM signal is given by 

2 i- r°o 

S F wm(u-,T,T u T 2 ,T 3 ) = - Im 


duj[ 

2 n 


r ‘°° rl / r°° rli / 

UUJ 2 j uoq e _j w r+jwiTi+^T2+j^T 3 


27T 


2 tt 


X S p (u>)£i(uJ 1 )S2((^2)^'3( u} 3)XfWm( ^1) W 2) ^ 3 ) ) (52) 


where the third order generalized susceptibility is 


1fU(- w . w I> w 2> w 3) = (53) 


where V = V + is the full matter transition operator. Note, that Eq. ( 1501) and Eq. f|52l) 
do not rely on the rotating wave approximation (RWA). 
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